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Abstract 

We examine Compton scattering of electrons on black body photons in the 
case where the electrons are highly relativistic, but the center of mass energy 
is small in comparison with the electron mass. We derive the partial lifetime of 
electrons in the LEP accelerator due to this form of scattering in the vacuum 
beam pipe and compare it with previous results. 
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I. INTRODUCTION 



Vacuum beam pipes of modern particle accelerators closely approach the ideal limit of 
a pipe completely devoid of gas molecules. However, even an ideal vacuum beam pipe in a 
laboratory at room temperature is filled with photons having an energy distribution given by 
Planck's law. Some time ago, Telnov0 noted that the scattering of electrons on these black 
body photons could be a significant mechanism for the depletion of the beam. This scattering 
of the electrons in the Large Electron Positron collider at CERN (LEP) on the black body 
radiation has been detected 2-4 . There is a long history of theoretical investigations on the 
scattering of high-energy electrons on black body photons, centering around the role this 
plays as an energy loss mechanism for cosmic rays, which is summarized by Blumenthal and 
Gouldl. More recently, Domenicoi and BurkhardtEl have considered this effect for the LEP 
experiments and the consequent limit on the beam lifetime by using numerical Monte Carlo 
methods. In view of the intrinsic interest of the problem of high-energy electron scattering on 
black body photons, we believe that it is worthwhile to present here a simplified calculation 
of the effect. We compute the total cross section analytically. The cross section as a function 
of the energy loss — which is the important quantity for the beam lifetime — is also done 
analytically except for a final straightforward numerical integration. Our calculations use 
relativistic invariant methods, and are thus of some pedagogical interest. 

Analytic computations can be performed because the problem involves two small dimen- 
sionless parameters. On the one hand, the electron of mass m has a very large laboratory 
energy E and it is ultrarelativistic, as characterized by the parameter m 2 /E 2 . (We use 
natural units in which the velocity of light c = 1, Planck's constant h — 1, and Boltzmann's 
constant k — 1, so that temperature is measured in energy units.) At LEP, m 2 /E 2 » 10 -10 . 
We shall neglect terms of order m 2 /E 2 . On the other hand, the temperature T of the black 
body radiation is very small in comparison with the electron mass m. Thus, although the 
electron is ultrarelativistic, the energy in the center of mass of the electron-photon system 
is still small in comparison with the electron mass. The head-on collision of a photon of 
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energy T with an relativistic electron of energy E produces, with the neglect of the electron 
mass, the squared center-of-mass energy 4ET. We shall use the dimensionless parameter 
(which gives an average value) 



At LEP, s « 1(T 2 . Thus it is a good approximation to use the non-relativistic limit in the 
center of mass, with the relativistic Compton cross section replaced by its constant, non- 
relativistic Thomson limit. To assess this approximation, we shall also compute the first 
corrections in s. 

In the next section we use simple relativistic techniques to compute the total cross section 
for the scattering of an ultrarelativistic electron on the Planck distribution of black body 
photons. The third section describes the more detailed calculation needed for the cross 
section in which the electron loses an energy greater than AE. If the energy loss AE in an 
electron-photon collision is too large, the electron's motion falls outside of the acceptance 
parameters of the machine. At LEP this happens when AE/E is greater than about 1%. 
As we shall see, this means that even if the beam were in a perfect vacuum, it would decay 
with a half life of about two days. The vacuum in the LEP accelerator is so good that the 
beam scattering of the black body photons is, in fact, the primary mechanism for beam loss 
when the machine is run with a single beam. Scattering on the residual gas in the beam 
pipe gives a considerably longer half life of about six days.H When the machine is run in the 
usual mode with two beams for e + e~ collision experiments, beam-beam collisions reduce 
the beam half life to about 14 hours.! 



In the general case of an electron scattering off some photon distribution, the scattering 
rate T in the electron's rest frame may be computed using the formula 



s = 2ET/m 2 . 



(1.1) 



II. TOTAL SCATTERING RATE 




(2.1) 
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where /(k) is the photon phase-space density [with the normalization defined such that j° 
in Eq. (|2.7|) is the photon number density] as a function of the photon momentum and cr(k) 
is the scattering cross section, which is similarly a function of the photon momentum. Here 
all quantities are evaluated in the electron's rest frame as indicated by the over bar. This 
scattering rate may be viewed as a time derivative 

— dn 

r = *- < 2 - 2) 

where r is the time in the electron's rest frame. Since numbers are Lorentz invariant and 
r may be defined to be the invariant proper time of the electron, the rate dn/dr is, in 
fact, a Lorentz invariant. Thus, if the integral on the right-hand side of the rate formula 
( [2.1D is written in a Lorentz invariant manner, we can immediately evaluate the rate in the 
laboratory frame. In the lab frame, the electron moves with four momenta 

= m — , (2.3) 
dr 

where z^ir) is the world line of this particle, its space-time position as a function of proper 
time, and m is the electron mass. This gives the familiar time-dilation formula 

tt = dz^ = p^ = E^ 
dr dr mm 1 

where E is the electron's total energy. Thus the scattering rate V in the lab frame may be 
easily evaluated using 

dn mdn 

~dt = ~E~drr ' (2 ' 5) 

In the non-relativistic limit, a may be replaced with the Thomson cross section, ot = 
8irr 2 /3, where r$ = e 2 /Anm is the classical electron radius. Since this is independent of k, 
the scattering rate may be rewritten as 

fo = <r T j°, (2.6) 

where 



is the photon number flux four vector. Since p^/m = (1,0) in the electron's rest frame, we 
may write this leading approximation, denoted with a subscript, as 

(*t) =-a T ]^, (2.8) 
\clt J m 

with the minus sign arising from our Lorentz metric convention in which the metric has 
the signature (— , +, +, +). The result ( |2.8|) is now in an invariant form which holds in any 
frame. With a thermal photon distribution in the lab frame, 

where to = k° is the photon energy, the photon number distribution is isotropic, and so only 
the number density component j° is nonvanishing. Thus, in the lab frame, 

r Q = ^a T f V - = a T f. (2.10) 
hi m 

The lab photon number density obtained from integrating Q2.7| ) with the distribution ( p.9| ) 
is the familiar result 

j0= 2C^) T 3 ; (2n) 
in which £(3) = 1.202 ... is the Riemann zeta function. Thus, T is given by 

r = ^TV- (2-12) 

The first order relativistic correction to this result is obtained with the use of the corrected 
cross section 

a = a T (l+ 2 -4). (2.13) 

Note that the product kp = k^p^ of the two four-momenta is negative with our metric. 
Because a is no longer independent of k, the corresponding form of Eq. ( |2.10| ) is 



m 



(d 3 k) —kp ( 2pk^ 



ri = 7^ / ^kT^f(k) 1 + 



E J (2vr) 3 k°m J v ' \ m? 
E \ m 2 



-^a T (f Pli + —T^ PlM p v ) , (2.14) 



where 



is the stress-energy tensor of the photons. Due to the isotropy of the thermal photons in the 
lab frame, T^ u has no off diagonal components, and it is also traceless because the photon 
is massless, k^k^ = 0. Therefore, in the lab frame, 

T^p,p u = (e 2 + l|p 2 |) T 00 = \e- [l - ^ T 00 . (2.16) 

The m? / E 2 term is very small, and it may be neglected. Integrating over the photon 
distribution in Eq. ( |2.15| ) gives the well-known black body energy density 

T oo = 6^4) T 4 ; (217) 
where £(4) = 7r 4 /90 = 1.082 .... This yields the corrected scattering rate 



1 1 = — — 1 a T 

7T^ 



C(4) 
1 -4s SV ; 



(2.18) 



C(3). 

For the temperature in the LEP beam pipe we take T = 291 K = 0.0251 eV, which is 
about room temperature. This gives the leading rate r = 9.98 x 10~ 6 s _1 corresponding to 
the mean life To = l/To = 28 hr. A typical LEP beam energy E = 46.1 GeV is just above 
half the Z° mass — within the width, but on the high side of resonance curve. Together with 
the previous value of the temperature, this gives s = 0.00886, and the first-order corrected 
rate Ti = 9.66 x 10~ 6 s _1 , which is about 3% smaller than the leading rate. This gives a 
mean life T\ = 1 /Y\ = 29 hr. 



III. RATE WITH ENERGY LOSS 



The calculation of the scattering rate in which the electron loses an energy greater than 
AE is facilitated by going back to the basic formula§ that expresses the rate in terms of 
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Lorentz invariant phase space integrals, an energy-momentum conserving 5 function, and 
the square of the scattering amplitude |T| 2 . The total electron scattering rate as observed 
in the lab frame reads 

r " h I {§> Iffih! W + <->- rtPT . (3-D 

where p and p' are the initial and final electron four momenta, k and k' the initial and final 
photon four momenta, with E = p° , E' = p , uj = k° , to' = k'° the time components of 
these four vectors. Except for the initial factor of 1/2E which is the lab energy of the initial 
electron and which converts the invariant proper time into the lab time, the right-hand side 
of this expression is a Lorentz invariant. The problem proves to be greatly simplified if the 
integrals are evaluated in the rest frame of the electron, because Compton scattering of a 
photon on an electron at rest has a very simple nonrelativistic limit. This complicates the 
initial photon distribution, but, if we introduce a four vector (3^, whose time component 
in the lab frame is one over the temperature of the photon distribution and whose spatial 
components are zero in the lab frame, the distribution in an arbitrary frame still has the 
simple form 

From the definition of —(3^/3^ = 1/T 2 and — = E/T, because multiplication by (3^ 
selects out the time component in the lab frame. In the electron rest frame, /3 M therefore 
takes on the value 

/^(^ (3.3) 

where E and p are taken in the lab frame. We have not yet taken into account the lower 
bound on the electron energy loss in the lab frame. Because multiplication by (5 11 selects the 
time component in the lab frame, this limit may be instituted by the inclusion of an "energy 
loss" step function in the integrand of Eq. (|3.1|) , 

AE\ 







(-/^(Pm-p;)-— ) , (3-4) 



where the 1/T in the second term has been included to compensate for the factor of 1/T 
which the first term has picked up by being multiplied by Using the identity 

Iwt^=Jwf e(E '- m)2wHp,2+m2) ' (3 - 5) 

the final electron four momentum p' may be integrated over to leave 

r < AB) = 21 / H2>> / - ^(-^ + - 2 "'> 

e(-/3fc'-(^-/*))|If. (3.6) 

We do the k' integral in spherical coordinates and take the z axis to be parallel to k, 
with 9 the angle between these two vectors. The angle 9 is the photon scattering angle in 
the electron rest frame, and 

-kk' = ujuj'(1- cos 9). (3.7) 

The 5 function can now be solved for uj' to yield 

5{-2muj + 2mw' - 2kk') = - i —5 (uj 1 ^ ] , (3.8) 

v ' 2(m + uj(l -cos 9)) \ m + oj(l-cos9)J v ' 

which requires that uj' < uj and thus makes the 9(uj — uj') step function redundant. The 
scattered photon energy uj' given by the 5 function is, of course, just the Compton energy. 
To deal with the "energy loss" step function, we note that the Lorentz transformation from 
the lab frame to the initial electron rest frame turns the lab frame isotropic black body 
photon distribution into a very narrow pencil in the electron rest frame in which we are now 
working. Thus the initial photon distribution is sharply peaked about the average value 

(k^/uj) = p*l (3° . (3.9) 

Hence we can approximate 

8° 

- f3k' ~ —{-kk') = P°uj'(1 - cos 9) . (3.10) 

UJ 

To verify this and assess the order of accuracy, we define the average more precisely by 
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( X ) = {2 ; LT . (3.ii) 

Then, by virtue of the relativistic invariance of this definition, 

(pr) = (V/r - ~/3V) MP 2 ) , (3-i2) 

since (3 X is the only four vector available and k^k^ = 0. Remembering that u = k°, this 
presents the squared fluctuation about the average as 

Wj' ( } 

where 

-/3 2 m 2 , . 

B = ~ . (3.14) 

Thus the deviations away from our approximation may be neglected because they involve 
the very small quantity m 2 /E 2 . Using this approximation for —k'f3 simplifies the "energy 
loss" step function to 

where the — (3k in the original step function has been neglected because it is much less than 
AE/T. Inserting the value of uo' given by the energy- conserving 6 function ( |3.<j| ) into the 
step function gives 

\m + uj(l-cos6)~ E l-cos0/ V ~ E' 1 - cos#J ' ^' ' 

where on the right hand side we have solved for u and defined 

E' = E — AE , (3.17) 

which is the maximum final electron energy in the lab frame. 

We perform the k integral in spherical coordinates, with the polar angle x taken to be 
the angle between k and /?, so that 
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- k(3 = oj((3°- |0| cos x)- (3.18) 
We rewrite the angular integral for k in terms of an integration over k/3 by noting the limits 

-kp<u(^+\0\)^2u;P o = ^, (3.19) 

and 

Thus, with the neglect of order m 2 /E 2 corrections, and remembering that \/3\ ~ E/mT, 

rl m T r 2ujE/mT 

/ dcos X = —^ d{-k(3). (3.21) 

J-l UOE Juim/2TE 

We shall do the —k/3 integral last, due to its dependence on the initial photon distribution. 
In order to interchange the order of the uj and — k/3 integrations, we note that the lower 
limit on —k/3, 

gives the upper bound on uj, 

uj < 2TE{ -~ k ^ = s (-k(3)m . (3.23) 
m 



The upper limit on —k/3, 

mT ' 



k/3<^i, (3.24) 



gives the lower bound on uu, 



( y —kj3)mT . . m 



2 



uj > = s (-k/3)m— . (3.25) 

In view of the extreme smallness of m 2 / E 2 , we may replace this lower limit by uj = 0. Hence, 
switching the order of integration gives 



/ duo d(-k/3) = d{-k(3) duo. (3.26) 

JO Ju>m/2TE JO Jo 
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We perform this reversal of integrals, do the two trivial azimuthal integrals, and do the uj' 
integral using the 5 function to obtain 

[m + w(l-cos9)]2 |T|2 6 k" ~ ~Wl-cme) ' (3 ' 2?) 
To work out the integrals which appear here, it is convenient to first introduce the 
appropriate, dimensionless variables, 

x = —k/3 , z — \ — cos#, v = — — ^— — , (3.28) 

s{—kp)m 

and define 

»=w- (329) 

With this new notation, we have 

T 3 r°° r 2 r 1 v f 2n\ 

W = **™l*l*Ji + ^W 9 ("-»)■ (3 ' 30) 

The final step function provides the lower limit v = 2u/xz. This lower limit must not exceed 
the upper limit v — 1. Hence we must have condition z > 2u/x on the z integration. But 
again, this must not exceed the upper limit z — 2. Thus x > u, and imposing all these limits 
gives 

T 3 r°° r 2 r 1 v 

r(A£) = - 9 . .„ / dxx 2 f{x) dz dv- ^l T ! 2 - ( 3 - 31 ) 

v 1 4m 2 (27r) 3 k v y kujx ku/xz (1 + svxz) 2 1 1 v 1 

This will be evaluated in the nonrelativistic limit, keeping first order corrections in s. 
The exact squared amplitude differs from its nonrelativistic limit 

|T| 2 = 2e 4 (l + cos 2 6) = 2e 4 (2 - 2z + z 2 ) (3.32) 

by corrections of order uouo' /m 2 . These corrections involve s 2 and are thus negligible. To 
first order in s, 

(i+Lz)* * 1 - 23 "**- (3 - 33) 
ii 



The z and v integrations are now straightforward. We express the result as 



F(AE) = r 



F (u) - As^-F 1 (u) 



(3.34) 



where To is the approximate total scattering rate from Eq. ( [2.12| ). The straightforward 
integrations give 



1 r°° dx ( /x\ 2u 3 \ 

= 203) /. —1 {* ~ 3UX + 3U ln U + -)' < 3 - 35 ) 



an 



d0 



oo 



V ; 6C(4) Ju e x -l\ 2 \u) 2x ) K ' 

It can be seen that in the AE limit, -F (0) = -Fi(O) = 1, so T(AE) reduces to the result 
( CT ) for Tl 

At this stage, one must resort to numerical calculations to evaluate the integrals. How- 
ever, analytic calculations of the energy weighted moments of the distribution can still 
be performed. The simplest of these is the average energy loss observed in the lab frame, 
(E — E'). Because r(A£') describes the rate due to all scattering events where E — E'> AE, 
this average value may be computed by 

where Ti is the total scattering rate including the first correction in s given by Eq. ( p,18|) . 
Changing variables to u and integrating by parts gives 



oo 



, , V(AE)dAE 
(E-E')= du K > , 3.38 
lo 1 1 du 



with, in view of Eq. ( |3.29| ) 



dAE _ 2sE 

~d^~ {l + 2suf (3 ' 39) 

Inserting the expressions for T(AE) and Ti into the integral and expanding in powers of s 
gives 
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(E - E') = 2s E J™ du |f («) - 4s 



u F (u)-^F (u) + ^F l{ a) 



+ 0(s 2 ) . (3.40) 



C(3) uv ' C(3)' 

Inserting the expressions for the F's from Eq. (|3.35| ) and Eq. (|3.36|) and interchanging the 
order of the x and u integrals, the integrals may be evaluated analytically, and we find that 

,C(4) 63C(5)" 



(E-E') = 2sE^h+s 



C(3) 5 C(4) 



+ 0(s' 



1.80 sE 1 -8.5 s 



(3.41) 



To check that no mistakes have been made in our calculation of Fq{u) and F\{u) given in 
Eq. (|3.35|) and Eq. Q3.36Q , we have independently evaluated the average energy loss (E — E') 
starting from Eq. ( |3.6|) and only making the small s approximation towards the end of the 
calculation. We find the same result with this different method. 

The integrals in the definitions ( |3.35|) and ( |3.36|) of the functions F${u) and Fi(u) have 
been calculated numerically, and the results are displayed in Fig. [L] and Fig. ^]. As a check 
on this numerical result, we have used it to evaluate the integrals in Eq. ( ^.40[ ) numerically, 
and the results agree with the analytic expression ( |3.41| ) to within 0.2 percent. 

We may compare our calculations with the those of Domenicoi who employed a Monte 
Carlo method. He used the values E = 46.1 GeV and T = 291 K (which we have previously 
employed) that give s = 0.0089. He also took AE = 0.012.E which places u = 0.69. 
Numerical integration gives F (0.69) = 0.44 and i<i(0.69) = 0.83, and from these values we 
calculate a mean beam lifetime of 64 hours to zeroth order in the nonrelativistic limit, and 
of 68 hours when the first order relativistic corrections are included. This is to be compared 
with Domenico's value of 90 hours for the same input parameters. We do not understand 
the reason for this discrepancy. 

We may also compare our results with those of H. Burkhardt@ who also used a Monte 
Carlo method. Burkhardt uses parameters that are slightly different than those used by 
Domenico, namely E = 45.6 GeV and T = 295 K. These parameters yield the same s = 
0.0089. However, since the overall rate scales as T 3 , which is 4% larger with Burkhardt's 
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FIG. 1. The dimensionless Fq(u) defined in Eq. (3.35) as a function of the dimensionless variable 
u defined in Eq. (3.29). 

temperature, our value of the beam lifetime for AE = 0.012.E with his parameters is reduced 
from 68 to 65 hours. Burkhardt finds 83 hours. (To compare with Domenico, we note that 
modifying his result of 90 hours by the 4% change in T 3 produces 86 hours.) Burkhardt 
and Kleissi also state that the average fractional energy loss (E — E')/E is 1.1% for this 
value of s, but our result ( 3.41]) gives the larger value 1.5% corresponding to our shorter 



beam lifetime. Again, we we can only state that we do not understand the reason for these 
discrepancies. 
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FIG. 2. The dimensionless Fi(u) defined in Eq. (3.36) as a function of u . 
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